TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway.

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

Genome-wide analysis of copy-number variation in humans with cleft lip and/or cleft palate identifies COBLL1, RIC1, and ARHGEF38 as clefting genes.

Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.

Hypersensitivity of Zebrafish htr2b Mutant Embryos to Sertraline Indicates a Role for Serotonin Signaling in Cardiac Development.

ABSTRACT: Selective serotonin reuptake inhibitors (SSRIs) are antidepressants prescribed in 10% of pregnancies in the United States. Maternal use of SSRIs has been linked to an elevated rate of congenital heart defects, but the exact mechanism of pathogenesis is unknown. Previously, we have shown a decrease in cardiomyocyte proliferation, left ventricle size, and reduced cardiac expression of the serotonin receptor 5-HT 2B in offspring of mice exposed to the SSRI sertraline during pregnancy, relative to offspring of untreated mice. These results suggest that disruption of serotonin signaling leads to heart defects. Supporting this conclusion, we show here that zebrafish embryos exposed to sertraline develop with a smaller ventricle, reduced cardiomyocyte number, and lower cardiac expression of htr2b relative to untreated embryos. Moreover, zebrafish embryos homozygous for a nonsense mutation of htr2b ( htr2bsa16649 ) were sensitized to sertraline treatment relative to wild-type embryos. Specifically, the ventricle area was reduced in the homozygous htr2b mutants treated with sertraline compared with wild-type embryos treated with sertraline and homozygous htr2b mutants treated with vehicle control. Whereas long-term effects on left ventricle shortening fraction and stroke volume were observed by echocardiography in adult mice exposed to sertraline in utero, echocardiograms of adult zebrafish exposed to sertraline as embryos were normal. These results implicate the 5-HT 2B receptor functions in heart development and suggest zebrafish are a relevant animal model that can be used to investigate the connection between maternal SSRI use and elevated risk of congenital heart defects.

TFAP2 paralogs facilitate chromatin access for MITF at pigmentation and cell proliferation genes.

In developing melanocytes and in melanoma cells, multiple paralogs of the Activating-enhancer-binding Protein 2 family of transcription factors (TFAP2) contribute to expression of genes encoding pigmentation regulators, but their interaction with Microphthalmia transcription factor (MITF), a master regulator of these cells, is unclear. Supporting the model that TFAP2 facilitates MITF's ability to activate expression of pigmentation genes, single-cell seq analysis of zebrafish embryos revealed that pigmentation genes are only expressed in the subset of mitfa-expressing cells that also express tfap2 paralogs. To test this model in SK-MEL-28 melanoma cells we deleted the two TFAP2 paralogs with highest expression, TFAP2A and TFAP2C, creating TFAP2 knockout (TFAP2-KO) cells. We then assessed gene expression, chromatin accessibility, binding of TFAP2A and of MITF, and the chromatin marks H3K27Ac and H3K27Me3 which are characteristic of active enhancers and silenced chromatin, respectively. Integrated analyses of these datasets indicate TFAP2 paralogs directly activate enhancers near genes enriched for roles in pigmentation and proliferation, and directly repress enhancers near genes enriched for roles in cell adhesion. Consistently, compared to WT cells, TFAP2-KO cells proliferate less and adhere to one another more. TFAP2 paralogs and MITF co-operatively activate a subset of enhancers, with the former necessary for MITF binding and chromatin accessibility. By contrast, TFAP2 paralogs and MITF do not appear to co-operatively inhibit enhancers. These studies reveal a mechanism by which TFAP2 profoundly influences the set of genes activated by MITF, and thereby the phenotype of pigment cells and melanoma cells.

Generating Zebrafish RNA-Less Mutant Alleles by Deleting Gene Promoters with CRISPR/Cas9.

Danio rerio (zebrafish), traditionally used in forward genetic screens, has in the last decade become a popular model for reverse genetic studies with the introduction of TALENS, zinc finger nucleases, and CRISPR/Cas9. Unexpectedly, homozygous frameshift mutations generated by these tools frequently result in phenotypes that are less penetrant than those seen in embryos injected with antisense morpholino oligonucleotides targeting the same gene. One explanation for the difference is that some frameshift mutations result in nonsense-mediated decay of the gene transcript, a process which can induce expression of homologous genes. This form of genetic compensation, called transcriptional adaptation, does not occur when the mutant allele results in no RNA transcripts being produced from the targeted gene. Such RNA-less mutants can be generated by deleting a gene's promoter using a pair of guide RNAs and Cas9 protein. Here, we present a protocol and use it to generate alleles of arhgap29b and slc41a1 that lack detectable zygotic transcription. In the case of the arhgap29b mutant, an emerging phenotype did not segregate with the promoter deletion mutation, highlighting the potential for off-target mutagenesis with these tools. In summary, this chapter describes a method to generate zebrafish mutants that avoid a form of genetic compensation that occurs in many frameshift mutants.

Vesicular monoamine transporter 2 (SLC18A2) regulates monoamine turnover and brain development in zebrafish.

AIM: We aimed at identifying potential roles of vesicular monoamine transporter 2, also known as Solute Carrier protein 18 A2 (SLC18A2) (hereafter, Vmat2), in brain monoamine regulation, their turnover, behaviour and brain development using a novel zebrafish model. METHODS: A zebrafish strain lacking functional Vmat2 was generated with the CRISPR/Cas9 system. Larval behaviour and heart rate were monitored. Monoamines and their metabolites were analysed with high-pressure liquid chromatography. Amine synthesising and degrading enzymes, and genes essential for brain development, were analysed with quantitative PCR, in situ hybridisation and immunocytochemistry. RESULTS: The 5-bp deletion in exon 3 caused an early frameshift and was lethal within 2 weeks post-fertilisation. Homozygous mutants (hereafter, mutants) displayed normal low locomotor activity during night-time but aberrant response to illumination changes. In mutants dopamine, noradrenaline, 5-hydroxytryptamine and histamine levels were reduced, whereas levels of dopamine and 5-hydroxytryptamine metabolites were increased, implying elevated monoamine turnover. Consistently, there were fewer histamine, 5-hydroxytryptamine and dopamine immunoreactive cells. Cellular dopamine immunostaining, in wild-type larvae more prominent in tyrosine hydroxylase 1 (Th1)-expressing than in Th2-expressing neurons, was absent in mutants. Despite reduced dopamine levels, mutants presented upregulated dopamine-synthesising enzymes. Further, in mutants the number of histidine decarboxylase-expressing neurons was increased, notch1a and pax2a were downregulated in brain proliferative zones. CONCLUSION: Lack of Vmat2 increases monoamine turnover and upregulates genes encoding amine-synthesising enzymes, including histidine decarboxylase. Notch1a and pax2a, genes implicated in stem cell development, are downregulated in mutants. The zebrafish vmat2 mutant strain may be a useful model to study how monoamine transport affects brain development and function, and for use in drug screening.

The opioid antagonist naltrexone decreases seizure-like activity in genetic and chemically induced epilepsy models.

OBJECTIVE: A significant number of epileptic patients fail to respond to available anticonvulsive medications. To find new anticonvulsive medications, we evaluated FDA-approved drugs not known to be anticonvulsants. Using zebrafish larvae as an initial model system, we found that the opioid antagonist naltrexone exhibited an anticonvulsant effect. We validated this effect in three other epilepsy models and present naltrexone as a promising anticonvulsive candidate. METHODS: Candidate anticonvulsant drugs, determined by our prior transcriptomics analysis of hippocampal tissue, were evaluated in a larval zebrafish model of human Dravet syndrome (scn1Lab mutants), in wild-type zebrafish larvae treated with the pro-convulsant drug pentylenetetrazole (PTZ), in wild-type C57bl/6J acute brain slices exposed to PTZ, and in wild-type mice treated with PTZ in vivo. Abnormal locomotion was determined behaviorally in zebrafish and mice and by field potential in neocortex layer IV/V and CA1 stratum pyramidale in the hippocampus. RESULTS: The opioid antagonist naltrexone decreased abnormal locomotion in the larval zebrafish model of human Dravet syndrome (scn1Lab mutants) and wild-type larvae treated with the pro-convulsant drug PTZ. Naltrexone also decreased seizure-like events in acute brain slices of wild-type mice, and the duration and number of seizures in adult mice injected with PTZ. SIGNIFICANCE: Our data reveal that naltrexone has anticonvulsive properties and is a candidate drug for seizure treatment.

Cooperation between melanoma cell states promotes metastasis through heterotypic cluster formation.

Melanomas can have multiple coexisting cell states, including proliferative (PRO) versus invasive (INV) subpopulations that represent a "go or grow" trade-off; however, how these populations interact is poorly understood. Using a combination of zebrafish modeling and analysis of patient samples, we show that INV and PRO cells form spatially structured heterotypic clusters and cooperate in the seeding of metastasis, maintaining cell state heterogeneity. INV cells adhere tightly to each other and form clusters with a rim of PRO cells. Intravital imaging demonstrated cooperation in which INV cells facilitate dissemination of less metastatic PRO cells. We identified the TFAP2 neural crest transcription factor as a master regulator of clustering and PRO/INV states. Isolation of clusters from patients with metastatic melanoma revealed a subset with heterotypic PRO-INV clusters. Our data suggest a framework for the co-existence of these two divergent cell populations, in which heterotypic clusters promote metastasis via cell-cell cooperation.

BRN2 is a non-canonical melanoma tumor-suppressor.

While the major drivers of melanoma initiation, including activation of NRAS/BRAF and loss of PTEN or CDKN2A, have been identified, the role of key transcription factors that impose altered transcriptional states in response to deregulated signaling is not well understood. The POU domain transcription factor BRN2 is a key regulator of melanoma invasion, yet its role in melanoma initiation remains unknown. Here, in a BrafV600E PtenF/+ context, we show that BRN2 haplo-insufficiency promotes melanoma initiation and metastasis. However, metastatic colonization is less efficient in the absence of Brn2. Mechanistically, BRN2 directly induces PTEN expression and in consequence represses PI3K signaling. Moreover, MITF, a BRN2 target, represses PTEN transcription. Collectively, our results suggest that on a PTEN heterozygous background somatic deletion of one BRN2 allele and temporal regulation of the other allele elicits melanoma initiation and progression.

MITF reprograms the extracellular matrix and focal adhesion in melanoma.

The microphthalmia-associated transcription factor (MITF) is a critical regulator of melanocyte development and differentiation. It also plays an important role in melanoma where it has been described as a molecular rheostat that, depending on activity levels, allows reversible switching between different cellular states. Here, we show that MITF directly represses the expression of genes associated with the extracellular matrix (ECM) and focal adhesion pathways in human melanoma cells as well as of regulators of epithelial-to-mesenchymal transition (EMT) such as CDH2, thus affecting cell morphology and cell-matrix interactions. Importantly, we show that these effects of MITF are reversible, as expected from the rheostat model. The number of focal adhesion points increased upon MITF knockdown, a feature observed in drug-resistant melanomas. Cells lacking MITF are similar to the cells of minimal residual disease observed in both human and zebrafish melanomas. Our results suggest that MITF plays a critical role as a repressor of gene expression and is actively involved in shaping the microenvironment of melanoma cells in a cell-autonomous manner.

Functional Characterization of a Novel IRF6 Frameshift Mutation From a Van Der Woude Syndrome Family.

BACKGROUND: Loss-of-function mutations in interferon regulatory factor-6 (IRF6) are responsible for about 70% of cases of Van Der Woude Syndrome (VWS), an autosomal dominant developmental disorder characterized by pits and/or sinuses of the lower lip and cleft lip, cleft palate, or both. METHODS: We collected a Chinese Han VWS pedigree, performed sequencing and screening for the causal gene mutant. Initially, species conservation analysis and homology protein modeling were used to predict the potential pathogenicity of mutations. To test whether a VWS family-derived mutant variant of IRF6 retained function, we carried out rescue assays in irf6 maternal-null mutant zebrafish embryos. To assess protein stability, we overexpressed reference and family-variants of IRF6 in vitro. RESULTS: We focused on a VWS family that includes a son with bilateral lip pits, uvula fissa and his father with bilateral cleft lip and palate. After sequencing and screening, a frameshift mutation of IRF6 was identified as the potential causal variant (NM.006147.3, c.1088-1091delTCTA; p.Ile363ArgfsTer33). The residues in this position are strongly conserved among species and homology modeling suggests the variant alters the protein structure. In irf6 maternal-null mutant zebrafish embryos the periderm differentiates abnormally and the embryos rupture and die during gastrulation. Injection of mRNA encoding the reference variant of human IRF6, but not of the frame-shift variant, rescued such embryos through gastrulation. Upon overexpression in HEK293FT cells, the IRF6 frame-shift mutant was relatively unstable and was preferentially targeted to the proteasome in comparison to the reference variant. CONCLUSION: In this VWS pedigree, a novel frameshift of IRF6 was identified as the likely causative gene variant. It is a lost function mutation which could not rescue abnormal periderm phenotype in irf6 maternal-null zebrafish and which causes the protein be unstable through proteasome-dependent degradation.

Analysis of zebrafish periderm enhancers facilitates identification of a regulatory variant near human KRT8/18.

Genome-wide association studies for non-syndromic orofacial clefting (OFC) have identified single nucleotide polymorphisms (SNPs) at loci where the presumed risk-relevant gene is expressed in oral periderm. The functional subsets of such SNPs are difficult to predict because the sequence underpinnings of periderm enhancers are unknown. We applied ATAC-seq to models of human palate periderm, including zebrafish periderm, mouse embryonic palate epithelia, and a human oral epithelium cell line, and to complementary mesenchymal cell types. We identified sets of enhancers specific to the epithelial cells and trained gapped-kmer support-vector-machine classifiers on these sets. We used the classifiers to predict the effects of 14 OFC-associated SNPs at 12q13 near KRT18. All the classifiers picked the same SNP as having the strongest effect, but the significance was highest with the classifier trained on zebrafish periderm. Reporter and deletion analyses support this SNP as lying within a periderm enhancer regulating KRT18/KRT8 expression.

Mutations in GDF11 and the extracellular antagonist, Follistatin, as a likely cause of Mendelian forms of orofacial clefting in humans.

Cleft lip with or without cleft palate (CL/P) is generally viewed as a complex trait with multiple genetic and environmental contributions. In 70% of cases, CL/P presents as an isolated feature and/or deemed nonsyndromic. In the remaining 30%, CL/P is associated with multisystem phenotypes or clinically recognizable syndromes, many with a monogenic basis. Here we report the identification, via exome sequencing, of likely pathogenic variants in two genes that encode interacting proteins previously only linked to orofacial clefting in mouse models. A variant in GDF11 (encoding growth differentiation factor 11), predicting a p.(Arg298Gln) substitution at the Furin protease cleavage site, was identified in one family that segregated with CL/P and both rib and vertebral hypersegmentation, mirroring that seen in Gdf11 knockout mice. In the second family in which CL/P was the only phenotype, a mutation in FST (encoding the GDF11 antagonist, Follistatin) was identified that is predicted to result in a p.(Cys56Tyr) substitution in the region that binds GDF11. Functional assays demonstrated a significant impact of the specific mutated amino acids on FST and GDF11 function and, together with embryonic expression data, provide strong evidence for the importance of GDF11 and Follistatin in the regulation of human orofacial development.

Drug repositioning in epilepsy reveals novel antiseizure candidates.

Objective: Epilepsy treatment falls short in ~30% of cases. A better understanding of epilepsy pathophysiology can guide rational drug development in this difficult to treat condition. We tested a low-cost, drug-repositioning strategy to identify candidate epilepsy drugs that are already FDA-approved and might be immediately tested in epilepsy patients who require new therapies. Methods: Biopsies of spiking and nonspiking hippocampal brain tissue from six patients with unilateral mesial temporal lobe epilepsy were analyzed by RNA-Seq. These profiles were correlated with transcriptomes from cell lines treated with FDA-approved drugs, identifying compounds which were tested for therapeutic efficacy in a zebrafish seizure assay. Results: In spiking versus nonspiking biopsies, RNA-Seq identified 689 differentially expressed genes, 148 of which were previously cited in articles mentioning seizures or epilepsy. Differentially expressed genes were highly enriched for protein-protein interactions and formed three clusters with associated GO-terms including myelination, protein ubiquitination, and neuronal migration. Among the 184 compounds, a zebrafish seizure model tested the therapeutic efficacy of doxycycline, metformin, nifedipine, and pyrantel tartrate, with metformin, nifedipine, and pyrantel tartrate all showing efficacy. Interpretation: This proof-of-principle analysis suggests our powerful, rapid, cost-effective approach can likely be applied to other hard-to-treat diseases.

Identification of Isthmin 1 as a Novel Clefting and Craniofacial Patterning Gene in Humans.

Orofacial clefts are one of the most common birth defects, affecting 1-2 per 1000 births, and have a complex etiology. High-resolution array-based comparative genomic hybridization has increased the ability to detect copy number variants (CNVs) that can be causative for complex diseases such as cleft lip and/or palate. Utilizing this technique on 97 nonsyndromic cleft lip and palate cases and 43 cases with cleft palate only, we identified a heterozygous deletion of Isthmin 1 in one affected case, as well as a deletion in a second case that removes putative 3' regulatory information. Isthmin 1 is a strong candidate for clefting, as it is expressed in orofacial structures derived from the first branchial arch and is also in the same "synexpression group" as fibroblast growth factor 8 and sprouty RTK signaling antagonist 1a and 2, all of which have been associated with clefting. CNVs affecting Isthmin 1 are exceedingly rare in control populations, and Isthmin 1 scores as a likely haploinsufficiency locus. Confirming its role in craniofacial development, knockdown or clustered randomly interspaced short palindromic repeats/Cas9-generated mutation of isthmin 1 in Xenopus laevis resulted in mild to severe craniofacial dysmorphologies, with several individuals presenting with median clefts. Moreover, knockdown of isthmin 1 produced decreased expression of LIM homeobox 8, itself a gene associated with clefting, in regions of the face that pattern the maxilla. Our study demonstrates a successful pipeline from CNV identification of a candidate gene to functional validation in a vertebrate model system, and reveals Isthmin 1 as both a new human clefting locus as well as a key craniofacial patterning gene.

Zebrafish models of orofacial clefts.

Zebrafish is a model organism that affords experimental advantages toward investigating the normal function of genes associated with congenital birth defects. Here we summarize zebrafish studies of genes implicated in orofacial cleft (OFC). The most common use of zebrafish in this context has been to explore the normal function an OFC-associated gene product in craniofacial morphogenesis by inhibiting expression of its zebrafish ortholog. The most frequently deployed method has been to inject embryos with antisense morpholino oligonucleotides targeting the desired transcript. However, improvements in targeted mutagenesis strategies have led to widespread adoption of CRISPR/Cas9 technology. A second application of zebrafish has been for functional assays of gene variants found in OFC patients; such in vivo assays are valuable because the success of in silico methods for testing allele severity has been mixed. Finally, zebrafish have been used to test the tissue specificity of enhancers that harbor single nucleotide polymorphisms associated with risk for OFC. We review examples of each of these approaches in the context of genes that are implicated in syndromic and non-syndromic OFC. Developmental Dynamics 246:897-914, 2017. © 2017 Wiley Periodicals, Inc.

A Gene Implicated in Activation of Retinoic Acid Receptor Targets Is a Novel Renal Agenesis Gene in Humans.

Renal agenesis (RA) is one of the more extreme examples of congenital anomalies of the kidney and urinary tract (CAKUT). Bilateral renal agenesis is almost invariably fatal at birth, and unilateral renal agenesis can lead to future health issues including end-stage renal disease. Genetic investigations have identified several gene variants that cause RA, including EYA1, LHX1, and WT1 However, whereas compound null mutations of genes encoding α and γ retinoic acid receptors (RARs) cause RA in mice, to date there have been no reports of variants in RAR genes causing RA in humans. In this study, we carried out whole exome sequence analysis of two families showing inheritance of an RA phenotype, and in both identified a single candidate gene, GREB1L Analysis of a zebrafish greb1l loss-of-function mutant revealed defects in the pronephric kidney just prior to death, and F0 CRISPR/Cas9 mutagenesis of Greb1l in the mouse revealed kidney agenesis phenotypes, implicating Greb1l in this disorder. GREB1L resides in a chromatin complex with RAR members, and our data implicate GREB1L as a coactivator for RARs. This study is the first to associate a component of the RAR pathway with renal agenesis in humans.

Beyond MITF: Multiple transcription factors directly regulate the cellular phenotype in melanocytes and melanoma.

MITF governs multiple steps in the development of melanocytes, including specification from neural crest, growth, survival, and terminal differentiation. In addition, the level of MITF activity determines the phenotype adopted by melanoma cells, whether invasive, proliferative, or differentiated. However, MITF does not act alone. Here, we review literature on the transcription factors that co-regulate MITF-dependent genes. ChIP-seq studies have indicated that the transcription factors SOX10, YY1, and TFAP2A co-occupy subsets of regulatory elements bound by MITF in melanocytes. Analyses at single loci also support roles for LEF1, RB1, IRF4, and PAX3 acting in combination with MITF, while sequence motif analyses suggest that additional transcription factors colocalize with MITF at many melanocyte-specific regulatory elements. However, the precise biochemical functions of each of these MITF collaborators and their contributions to gene expression remain to be elucidated. Analogous to the transcriptional networks in morphogen-patterned tissues during embryogenesis, we anticipate that the level of MITF activity is controlled not only by the concentration of activated MITF, but also by additional transcription factors that either quantitatively or qualitatively influence the expression of MITF-target genes.

Identification of common non-coding variants at 1p22 that are functional for non-syndromic orofacial clefting.

Genome-wide association studies (GWAS) do not distinguish between single nucleotide polymorphisms (SNPs) that are causal and those that are merely in linkage-disequilibrium with causal mutations. Here we describe a versatile, functional pipeline and apply it to SNPs at 1p22, a locus identified in several GWAS for non-syndromic cleft lip with or without cleft palate (NS CL/P). First we amplified DNA elements containing the ten most-highly risk-associated SNPs and tested their enhancer activity in vitro, identifying three SNPs with allele-dependent effects on such activity. We then used in vivo reporter assays to test the tissue-specificity of these enhancers, chromatin configuration capture to test enhancer-promoter interactions, and genome editing in vitro to show allele-specific effects on ARHGAP29 expression and cell migration. Our results further indicate that two SNPs affect binding of CL/P-associated transcription factors, and one affects chromatin configuration. These results translate risk into potential mechanisms of pathogenesis.

TFAP2 paralogs regulate melanocyte differentiation in parallel with MITF.

Mutations in the gene encoding transcription factor TFAP2A result in pigmentation anomalies in model organisms and premature hair graying in humans. However, the pleiotropic functions of TFAP2A and its redundantly-acting paralogs have made the precise contribution of TFAP2-type activity to melanocyte differentiation unclear. Defining this contribution may help to explain why TFAP2A expression is reduced in advanced-stage melanoma compared to benign nevi. To identify genes with TFAP2A-dependent expression in melanocytes, we profile zebrafish tissue and mouse melanocytes deficient in Tfap2a, and find that expression of a small subset of genes underlying pigmentation phenotypes is TFAP2A-dependent, including Dct, Mc1r, Mlph, and Pmel. We then conduct TFAP2A ChIP-seq in mouse and human melanocytes and find that a much larger subset of pigmentation genes is associated with active regulatory elements bound by TFAP2A. These elements are also frequently bound by MITF, which is considered the "master regulator" of melanocyte development. For example, the promoter of TRPM1 is bound by both TFAP2A and MITF, and we show that the activity of a minimal TRPM1 promoter is lost upon deletion of the TFAP2A binding sites. However, the expression of Trpm1 is not TFAP2A-dependent, implying that additional TFAP2 paralogs function redundantly to drive melanocyte differentiation, which is consistent with previous results from zebrafish. Paralogs Tfap2a and Tfap2b are both expressed in mouse melanocytes, and we show that mouse embryos with Wnt1-Cre-mediated deletion of Tfap2a and Tfap2b in the neural crest almost completely lack melanocytes but retain neural crest-derived sensory ganglia. These results suggest that TFAP2 paralogs, like MITF, are also necessary for induction of the melanocyte lineage. Finally, we observe a genetic interaction between tfap2a and mitfa in zebrafish, but find that artificially elevating expression of tfap2a does not increase levels of melanin in mitfa hypomorphic or loss-of-function mutants. Collectively, these results show that TFAP2 paralogs, operating alongside lineage-specific transcription factors such as MITF, directly regulate effectors of terminal differentiation in melanocytes. In addition, they suggest that TFAP2A activity, like MITF activity, has the potential to modulate the phenotype of melanoma cells.